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! Abstract 

5h , Two recent test failures of Hypersonic Technology Vehicle 2 impose a strike to the increasingly growing enthusiasm, not only on 

the United States side. It is important to find out the exact failure reason, otherwise a solution is impossible. In this Note, we propose 
' a potential failure reason from the perspective of lateral stability analysis. We argue that the time variant pressure fluctuations, which 
are normally omitted in classical aircraft dynamics analysis, could not be neglected in dynamic analysis of hypersonic vehicles. To 
demonstrate the idea, a hypersonic model is imagined in this work and its aerodynamic parameters are estimated using fundamental fluid 
' principles. Pressure fluctuations are thereafter estimated by an empirical formula. A lateral dynamic equation is set up, taking those time 
I I variant fluctuations into account. The resulted equation is a Mathieu differential equation. Numerical solutions of this equation show that 
the inclusion of fluctuation terms generates more complicated dynamics and should be considered in flight controller design. 

^ Introduction 

^ypersonic vehicle has become a popular topic once again in aerospace community [1-5]. However, two recent test failures of Hyper- 
(_jonic Technology Vehicle (HTV-2) [6] impose a strike to the increasingly growing enthusiasm, not only on the United States side. It is 

' ^nportant to find out the exact failure reason, otherwise a solution is impossible. In this work, we propose a potential failure reason from 
^^e perspective of lateral stability. We argue that the time variant pressure fluctuations, which are normally omitted in classical aircraft 

^y namics analysis, could generate trouble for hypersonic vehicle dynamics. 

, Flight dynamics and control has been recognized as an important research topic since Wright brothers first flied to sky in 1903. 

Lateral and longitudinal stability are two classical flight dynamics and control problem. In the lateral stability problem, coupling 
^—dynamics was an interesting phenomenon that can produce either negative or beneficial stability for various aircraft design. According 

the released information, [6] serious dynamic coupling between yaw and roll causes, at least, the first test failure of HTV-2. Similar 
^l^oupling problems have been discovered in X-series aircrafts and space shuttle. A good summary of those NASA practices can be found 

the literature [7]. 

(SI . For the HTV-2 case, some one may argue that a coupling between pitch and lateral dynamics could exist. In contrast, we hypothesize 
.,_jfliat the failure could still happen in a classical way, that is, the dynamic divergence develops in lateral coordinates. However, we should 
(^t^ke new physical factors into account. It is well known that hypersonic flight generates huge acoustic pressure and suffers vibroacoustic 
(problems. It is interesting to know if the acoustic wave, or more generally, the fluctuating pressure, could impede the dynamic stability, 
methodology is introduced in this paper to analyze the effect of surface pressure fluctuations on lateral dynamic stability, by solving 
so-called Mathieu differential equation. It is worthwhile to point out that the proposed method should be generic, although the initial 
Jea comes from analyzing the HTV-2 case. 



^ Vehicle and Aerodynamic Information 

Pressure fluctuations in turbulent boundary layer [8] were mainly studied for vibroacoustics. Experimental investigations have been 
conducted at subsonic speed to measure surface pressure fluctuations on a smooth wall, which could induce panel fluttering and structural 
vibrations. 

Very limited experimental results have been documented for surface pressure fluctuations at hypersonic speed. Only some empirical 
equations have been proposed, based on preliminary experiments and analysis of a straight flat plate at Mach numbers 2, 5 and 8 
[9-11], where compressibility, heat transfer between mediums and viscous effects are incorporated. The equations are quite helpful. 
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summarizing physical insights and setting a good start point for dynamic analysis. For readers convenience, one equation adopted in 
this work is listed below. It is 

^™^" 1 + 0.13MI' ^' 

where Prms is the root mean square of pressure fluctuation on a flat plate, peV^ /2 is local dynamic pressure and is local Mach 
number. The development of the equations was based on experimental data and fluid principles. Details can be found in literatures [12]. 
The above formulation has been justified in various works [12]. It is worthwhile to emphasize that the formulation was presented for 
hypersonic flat plates at small angle of attack. A practical hypersonic vehicle (see Fig. 1) has more complicated shape and generally 
operates at a large angle of attack. The above formulation is adopted to assist preliminary analysis. 

It can be seen that surface pressure fluctuations can be estimated using Eq. (1) based on the knowledge of local flow field, which 
can be calculated given the information of aerodynamic profiles. However, the detailed specifications of HTV-2 are not released to 
the public. Hence, the following analysis is performed on our own model. Figure 1 shows the vehicle, which is drawn according to a 
beautiful artistic illustration of HTV-2. It can be seen that the cone shapes of HTV-2 are simplified to four flat surfaces (ABC, ABD, 
ACD, BCD in Fig. 2) in our model, which ease the following calculations. In addition, four rudders are deployed underneath our model. 
Two of these rudders control longitudinal dynamics and the other two can control lateral dynamics. This setup decouples roll and pitch 
actuators and thus simplify the ongoing controller design. The imaginary geometry specifications of our hypersonic vehicle are given in 
Fig. 2. The onboard reaction system is absent on this vehicle. 
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Figure 2: Geometrical specifications of our hypersonic model. 
For simplicity, the mass distribution of the vehicle model is assumed uniform and the density is 2000kg/m'^. The moment of inertia 

is 
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Theoretical methods are used to approximate local flowfield of our vehicle. A code in MATLAB is developed to automatically 
calculate local flowfield based on oblique-shock theory and Prandtl-Meyer theory, given the geometrical shapes, angle of attack (a), 
angle of sideslip {/3) and atmospheric conditions. Tables 1 and 2 give approximations of aerodynamic parameters at two operating 
conditions. The aerodynamic coefficients are nondimensionalized to freestream dynamic pressure. The rolling momentum L and the 
yaw moment N are also calculated based on the given geometrical specifications. 
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Table 1: Aerodynamic coefficients at height 20km, A/qo 10 and a Odeg. 



/3(deg) 




Cy 




Ci 


Cm 


Cn 


L(N-m) 


N(N-m) 
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-0.0016 
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2574.8780 
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-0.0015 
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Table 2: Aerodynamic coefficients at height 50km, Mrx 10 and a Odeg. 



/3(deg) 


Cx 


Cy 


c. 


Ci 


Cm 


Cn 


L(N-m) 


N(N-m) 
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0.0004 


14.2435 
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Figure 3: Relation between Prms ^nd q on (a) plane ABD and (b) plane ACD. 



Control and stability derivatives can be estimated by perturbing attitudes and rudder setups. For example, in the case of Table 
1, roll and yaw moment coefficients with respect to P is Lp ~ 513. 2N ■ m/dcg and Nj^ = 469. 6N • m/dcg. The coefficients are 
= 7.12N ■ m/deg and Nf} 6.51N • m/deg for the case of Table 2. 

Pressure fluctuations on the two lateral surfaces (ABD and ACD) can be calculated using the empirical formula (Eq. (3)) based 
on the fluid approximations. Figure 3 shows the results, where q is dynamic pressure of the freestream. It can be seen that pressure 
fluctuations rapidly diminish at high altitude. In contrast, the changes of pressure fluctuations are not so sensitive to q. In addition, 
Prms/q is identical at various heights according to the present formula. In this approximated investigation, the peak Pmis is about 
140dB, achieved at A/ = 5 and 10km height. It was reported in the literature [13] that the value of pressure fluctuations on the forward 
fuselage could reach 160dB, suggesting that the formula adopted here underestimates pressure fluctuations. However, the following 
analysis shows that even these potentially underestimated fluctuations could cause lateral dynamic instabilities. 



3 Lateral Stability Analysis 

Five-degree-of-freedom equations of motion are defined relative to body coordinates [14], 

U- Ixz Ixz Al Al Al AP AP 

P = i—j )qr+—r+—pq+—Cipl3 + —CisJa + —CisJr + —-f-Ci^P+—^ 

J^x J^x ^x J^x ^x ix ^Vlx 

1 = ( — -J )p-r+—r^-—p^ 2 — r + —Cm^Aa+ —C^,Je + —--C^^g, (3) 

y y y y y y zv ly 

,Ix-Iy. Ixz . Ixz , Ix^^^e Al „ Al Al AP AP 

r = ( p-^)Pg+ -^P- ^9''+ — S q+ C„i,l3+—ClsJr + —Cl6ja + —--Cn^r+-—CnpP, 

iz ^z ^z ^z ^z ^z ^z i Z '^Vlz 



— sin<^icos^? — rcosa + psina H 

V ^ mV 



3 A. a , ■ ^ 1^ 
a = — costpcosy + q — psma a, 

V mV ° 

where A = pV^S/2, S is aerodynamic surface area, c is the character length, V is speed of freestream, (p, r) are and angular rates, 
((/), 9) are Euler angles, I is wing span, 5a is rudder operating set, (C;, Cm, Gn) are roll, pitch, yaw momentum coefficients and m is 
mass. 

The following assumptions can be accepted to simplify the above equations: 

(1) a is time invariant during lateral maneuver, a = a^, 

(2) sina « a and cosa « 1; « 0; 

(3) terms divided by the relatively very large freestream speed {V) are neglected; 

(4) Ixz can be omitted. 

The simplified lateral-directional equation is 

p = Lpp/h ,r = Npp/I,,p^~r + aop. (4) 

It is easy to see 

^ = -r + aaP={-^+ao^)l3. (5) 
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The lateral dynamics can be examined by analyzing this equation. For example, inherent oscillation frequency (iOd) is determined as 



It should be noted that the aerodynamic coefficients in Eq. (5) actually consist of stationary and fluctuation parts, that is 

/3 = (-f^+ao^)/3 + (-^ + ao^)sinM)/3, (6) 

where the subscript (0) denotes stationary parts, and the subscript (/) denotes time variant parts whose frequency is oj/. Without loss 
of generality, only one single frequency is considered here. The time variant parts can be mainly from the surface pressure fluctuations, 
which are omitted in classical dynamic analysis for traditional flights. In this work, the above mentioned calculations of our hypersonic 
model suggest that the stationary aerodynamic coefficients are about one order of magnitude larger than the fluctuating pressure parts 
(e.g. see Table 1, Table 2 and Figure 3). As a result, we argue that it could be improper to neglect time variant terms in dynamic 
equations for hypersonic vehicles. We can take this argument one step further and hypothesize that a controller synthesized based on a 
dynamic model with the absence of those fluctuating coefficients may lead to HTV-2 test failures. 
To show the idea, the above equation is written in a succinct form 

/3 = -[a + 6sin(a;/t)]/3, (7) 

which is the so-called Mathieu differential equation [15], where a = —f^ — ao-f ™d b ~ —j^ ~ ao-p-. A MATLAB code is 
programmed to find the numerical solutions of Eq. (7). Given initial values of a, b, ojf, /?|t=o and $\t=o, the numerical solutions of 
P\t=T and $\t=T can be calculated using Runge-Kutta scheme. The stability of the dynamic equation can be examined by checking the 

value of A: = log]^o(/5|t=T /(3\t=o)/T, where T should be larger than the period of fluctuations (empirically, one order of magnitude). A 
positive value of k implies instability and vice versa. Contours of k define the stability region. 

Figure 4 shows the stability regions at Uf =lHz, lOHz and lOOHz, respectively. It can be seen that a vehicle can be unstable even if 
6 = 0, which is the case with the absence of fluctuating terms. It is also obvious, as negative values of a (and 6 = 0) lead the divergence 
of Eq. (7). The instability can be controlled with the introduction of external actuation. A careful vehicle design can also resolve issue 
by generating a desired a. We can find that exactly the same procedure has been performed after the first test failure of HTV-2. It was 
said that engineers adjusted center of gravity of the HTV-2 and diminished its angle of attack in the second flight trajectory (to adjust 
a). [6]. In addition, onboard reaction control inputs were enforced (to remedy a). Hence, the second test failure of HTV-2 should due 
to other causes. Figures 4(a)-(b) show one potential factor It can be seen that the stability of the system becomes extremely sensitive 
to the value of b, if the vehicle is suffered IHz pressure fluctuations, which can be caused by flow-induced noise and vibrations during 
hypersonic flight. Once the lateral dynamics becomes unstable, exactly the same phenomenon of HTV-2 flight anomaly can happen: a 
higher-than-predicted yawing motion presents, leading to the coupled rolling motion that finally fails the flight . [6] 

Figure 4 also show some interesting facts. Fluctuations at higher frequencies (lOHz and lOOHz) have little influence on the lateral 
stability, which counteract our intuition to some extent. It seems a good news as it is well known that the frequency of acoustic 
perturbation is beyond 20Hz and the frequency of vibration should be larger than 5Hz. However, Eq. (7) is a nonlinear equation (and so 
is hypersonic vehicle itself). High frequency parts can possibly affect low frequency parts as well. The dynamics of Eq. (7) can become 
intractably complicated if time variant terms of multi-frequencies are included. The further analysis is not conducted as we lack the 
knowledge of the spectrum of pressure fluctuations on hypersonic vehicles. 

4 Summary 

This paper proposes that pressure fluctuations could lead to lateral instability of hypersonic vehicles. A conceptual design imitating 
HTV-2 to some extent was produced. The related local flowfield was approximated using oblique-shock wave and Prandtl-Meyer 
theories. Pressure fluctuations over the vehicle surfaces were derived based on a previous empirical formulation for a hypersonic flat 
plate. It will be definitely helpful if either a more accurate and expensive computational fluid dynamics can be performed, or an empirical 
formulation for a hypersonic wedge is available. In addition, this paper only proposes one potential reason but without a solution. All 
those constitute ongoing works. It should be admitted that the initial idea develops from the analysis of HTV-2 anomaly. However, the 
analytical method developed in this Note should be generic. 
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A Notation 

The following symbols are used in this paper: 



A = pFV2,Pa 

c = the character length, m 

{Ci , Cm ,Cn)= roll, pitch, yaw momentum coefficients 

/ = moment of inertia, kg • 

I = the wing span, m 

L = roll moment, N • m 

M = Mach number 

m = mass of aircraft, kg 

TV = yaw moment, N • m 

{p,q,r) = roll, pitch and yaw rate, rad/s 

Prms = rms fluctuating pressure. Pa 

q = dynamic pressure. Pa 

S = aerodynamic surface area, 

{x, y, z) = coordinates, m 

V = speed of freestream, m/s 

a = angle of attack, deg 

(3 = angle of sideslip, deg 

{(j), 9) = Euler angles, deg 
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p = density of air, kg/ 
Subscripts 

e = properties local to aerodynamic shape 

/ = time variant part 

oo = freestream conditions 

= stationary parts 
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